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Abstract. Dense pockets of cold, molecular gas precede the formation of 
stars. During their infancy and later phases of evolution, stars inject con- 
siderable energy into the interstellar medium by driving shocks either due 
to ionising radiation or powerful winds. Interstellar shock-waves sweep 
up dense shells of gas that usually propagate at supersonic velocities. It 
is proposed, in this paper, to examine the possibility of dense structure- 
formation and perhaps, future protostar-formation, in a molecular cloud 
shocked by such a shell. Here I shall discuss results of a self-gravitating, 
3 -dimensional, high-resolution simulation using the smoothed particle hy- 
drodynamics. 
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1. Introduction 

Detailed mapping of nearby star-forming clouds at submm wavelengths within legacy 
surveys such as the SCUBA2 Gould-belt survey, the Spitzer, and more recent Her- 
schel survey of these clouds has revealed their richness in structure. Star-forming 
clouds occurring at different stages of evolution have various sizes, and shapes. Inter- 
stellar shocks profoundly affect gas dynamics on a large scale, and are likely to be 
crucial in producing the observed dense filamentary clouds. I propose to examine this 
hypothesis numerically. 

2. Description of the problem 

I shall consider a simple test case where a molecular cloud is shocked by an incident 
shell moving with a supersonic velocity. The cartoon in Figure 1 demonstrates the 
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I I putational domain. Shown here is a periodic 

box enclosing the cloud confined by an inter- 
Figure 1. A cartoon showing physical details of cloud medium, and a slab approaching it with a 
the slab-cloud system. See text for description. precoUision velocity. Vs. 



relevant physical details of the problem. The computational domain may be conveni- 
ently divided into three regions viz. the intercloud medium (ICM), the standing point 
at the cloud surface, and the interiors of the shocked cloud marked 1, 2 and 3, respect- 
ively. The slab, as shown in this cartoon, moves through the ICM from the left to the 
right with a precollision velocity, V, . It forms a standing shock, labelled 2 in Fig. 1, 
after colUding with the cloud surface. The impact of this collision generates reflected 
waves in region 1, while some straddle the cloud surface and others propagate within 
the post-collision cloud, the so called transmitted-shock, with a velocity, Vc, in re- 
gion 3. The reflected shock- waves have been marked with curly arrows in Figure 1. 
Pressure in each of the three regions will be denoted by p with a subscript 1, 2 or 3, 
commensurate with the notation introduced in Fig. 1 . 

Any discussion of this problem would be incomplete without appropriate refer- 
ence to the shock-dynamics, however, for want of space, I shall restrict myself to only 
quoting some important results which can be derived using the Rankine-Hugoniot 
jump conditions (e.g. Courant & Friedrichs 1956). Following the notation introduced 
above, the excess pressure due to a highly supersonic shock is 

^1Z£2. 2+1, where ;.2 = 21Z1; (i) 

Pi - po /^^ r + 1 

which for an adiabatic gas constant, y = 5/3, produces a pressure-excess of ~ 6. 
Pressure within the cloud due to the transmitted shock, p^, can be shown to be related 
to its preshock value, pq, as 



(1 -2M2)- 1/y 



(l-2M2)-l/y 



(2) 



where M+ is the Mach number for the incident shock. Observe that, relative to the in- 
cident Mach number, there is a fall in pressure at the surface of incidence, the standing 
point. Finally, the density within the shocked cloud compares with that in the external 
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medium as, 



P3 ^ (y- 1) -P3/po(3y- 1) 
po pa/poCr- i)-(3r- 1)' 



(3) 



For a finite value of this ratio, the velocity of transmitted shock, Vc, within the cloud 
can be calculated using the equation of continuity. The resulting effect will be dis- 
cussed in Section 5 below. 



We have used the Lagrangian, particle-based scheme. Smoothed particle hydrodynam- 
ics (Monaghan 1992), hereafter referred to as SPH. An SPH particle, in the strictest 
sense, is not a point particle but has a finite spatial extent defined by a quantity called 
the smoothing length, h. Each SPH particle, characterised by state properties of dens- 
ity, velocity and temperature, interacts with other particles through numerical viscos- 
ity. The density of a particle, p,, is calculated by summing over contributions from 
nearest neighbours, N„eibs, of that particle within a radius 2h, 



where r is the position vector of the particle with identifier, i. The density of the fluid, 
Ps, is then simply 



The search for the nearest neighbours, and the calculation of net force on an individual 

particle is done using a tree-algorithm (Barnes & Hut 1986). Local inhomogeneities 
are smoothed out using a kernel that has been normalised to unity. Gravity for van- 
ishingly small inter-particle separations is smoothed out using a modified kernel that 
introduces an infinitesimally small repulsive force between closely spaced particles. 
The Thomas-Couchman kernel, a special type of the M-4 cubic spline, ensures a fi- 
nite repulsive force between closely interacting particle-pairs (Thomas & Couchman 
1992). The numerical algorithm employed here, Seren, is a weU tested code (Hubber 
et al. 201 1), and includes all the features described above. 

The thermodynamic details of the problem are modelled using a simple barotropic 
equation of state (EOS), defined by Equation (6) below, that mimics the post-collision 
temperature-jump, Tp^, and cools off to the precoUision temperature, 7,^, at relatively 



3. Numerical scheme 



Neibs 




(4) 




(5) 
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higher densities. 





; lO-^^g cm-3 < p < 10-^2g cm-3 



- = fe/m) X 
P 



< p < 5 X 10"^^g cm"3 
(^);5xl0-22gcm-3 <p 

< 10"i^g 



(6) 



\ K ^\ io-i5g cm-' 
p > 10"^^g cm"^; 





where the adiabatic gas constant, y = 5/3, kg and m are respectively the Boltzmann 
constant, and mean molecular mass. 



We propose to test the above theory using the scheme outlined in Figure 2. Ordinary 
SPH particles, representing gas within the cloud and the slab, interact with each other 
via gravitational and hydrodynamic forces. The ICM confining the cloud, however, 
is represented by special particles that exert only thermal pressure on other particles. 
The entire assembly, including the slab, is enclosed in a self-wrapping periodic box, 
where the periodicity is limited to merely ghosting particles; in other words, a particle 
leaving through one face of the box enters from the opposite face. 

Elsewhere in the literature, we have discussed an ensemble of cases for different 
choices of resolution (Anathpindika & Bhatt 2011). In this paper though, the em- 
phasis being on the formation of dense structure within the shocked cloud, I will only 
discuss the simulation with highest resolution. The precollision cloud was modelled 
as a sphere having uniform density, and the probability distribution function (PDF) of 
particles in the virgin cloud peaking at the predicted density, ptrue, is shown in Figure 
3, which not only demonstrates the stabiUty of the precolUson cloud, but also the ab- 
sence of any spurious dense pockets. 
Important parameters Mass of the cloud, Mdd = 500 Mq; radius, Rcid= 0.5 pc; 
uniform temperature, Tiso = 15 K. The cloud, ICM, and the slab were represented by 
SPH particles having three different choices of mass, m. If /, j, and k are the respective 
identifiers of particles in each of the three regions then, (m,, m^, ot^)/[M0] = (2.2 X 
10"^, 2.53 X 10"^, 1.26 X 10"^), and the number of particles in each of these regions is 
respectively (2.4, 2.3, 7.3)xl0^ so that total number of particles, Nm = 1.2 x 10*. 
Resolution The smallest resolvable spatial scale in an SPH calculation, as noted 
above, is the average smoothing length, havg, so that the volume of a typical SPH 
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Figure 3. A histogram of the gas density within 
the precolhsion cloud. The relatively small 
spread around the expected density, p,rue, indic- 
ates the stability of this cloud. 



Figure 4. A plot showing the radial variation of 
pressure within the shocked cloud at different 
epochs. 




0.5 1 0.5 



0.5 1 



Figure 5. A rendered density plot showing a time sequence of the mid-plane of the shocked 
cloud, as seen along a direction orthogonal to the plane of the figure. Turbulence induced by 
the transmitted shock generates clumps and filaments visible in these plots. 



particle, Vsph 



2 1 Nm I 



?3 . 



(7) 



is the typical smoothing length of a particle in the test cloud. In the present case, 
havg ~ 1-5 X 10"^ pc, and IK = Ajjhaug ~ 13, where Aj = (na-^^/Gpy^^, is the length 
of the fastest growing unstable mode in a gas body, the Jeans length. The quantity K 
defines the number of SPH particles available to resolve the unstable mode, which in 
this case is ~ 6, and therefore satisfies the Truelove criterion of spatially resolving the 
instability (Truelove et al. 1998). 



5. Result 



The post-collision reflected shock-wave as well as the transmitted shock within the 
cloud is evident from the plot in Figure 4, that shows the radial variation of pressure 
at different epochs. The incidence of the front surface of the slab on the cloud causes 
a jump in pressure at its surface, defined by Equation (1), followed shortly by another. 
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relatively smaller, jump due to the rear surface of the slab shocking the cloud (red- 
curve). The green curve shows the pressure distribution at a slightly later epoch when 
the pressure -jump moved leftward relative to its position in the red curve, suggesting 
an outward propagation of a wave, i.e. the reflected shock moving in the ICM. It 
can be identified as the flared region around the shocked cloud shown in the rendered 
density plot of Figure 5. 

The weaker transmitted shock that is most important in the evolution of the post- 
collision cloud, initially causes slight contraction of the cloud, evident from the col- 
lection of denser gas in a ring on its periphery. Relative to the incident shock, gas 
within the cloud is at a lower pressure manifested by a trough in the pressure distribu- 
tion at the surface of incidence, defined by Equation (2). Inward propagation of this 
shock generates a turbulent velocity field within the cloud that soon produces fractal 
structure, in other words the appearance of relatively dense pockets of gas interspersed 
with rarefied regions, the so called holes. Structure within the shocked cloud grows 
on a timescale comparable to the crossing time, tc ~ IRddlVs, of the precoUision slab 
which is much shorter than the growth time, tg, of the fastest growing unstable mode 
of length, A ~ (Trv^^j-IGpcidy^^, where Vgfj = ao + ''e; Vc ~ Vs/M+ ~ 0.4 km/s, and 
ao(r = 2QK) ~ 0.27 km/s, so fliat A -0.1 pc, which imphes, tg ~ 0.14 Myr. 

The simulation was terminated when the slab, having traversed the width of the 
cloud, reached the opposite face of the periodic box. Consequently, further evolution 
of clumps and other contiguous filaments could not be investigated in this work. We 
have observed that the incident slab, after colliding with the cloud, sufl'ers severe 
ablation by the time it reaches the other end of the cloud. The transmitted shock, 
relative to the incident slab, propagates at a much lower velocity within the post- 
collision cloud, evident from the rendered density plots in Figure (5). However, as 
argued in Section 2 above. Equations (8) below show that the transmitted shock in a 
cold, dense gas could still be supersonic, and therefore, will likely have a significant 
effect on the internal structure of the cloud. The pressure behind the slab changes on 
a very short timescale, the dynamical timescale, td, of the slab which can be shown to 
be 

where 

and M = 'VdaQ, is the Mach number for the transmitted shock; see paper I for the 
derivation. In the present case, td ~ iO''^ Myrs « tc, implying, the dynamical prop- 
erties behind the shock must indeed change rapidly. In this Ught the problem under 
consideration here, in effect, reduces to one of a weak-shock impinging on a cloud or 
conversely, a shock interacting with a large cloud. 
The probability distribution function (PDF) The CMF has often been suggested 
as the likely progenitor of the steUar IMF, and the likely relation between the two dis- 
tributions has been exannined by several authors (e.g. Anathpindika 2011, HatcheU & 
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Figure 6. The density PDF for gas in the 
shocked cloud after the slab has traversed the 
cloud. Bimodal nature of the density distribu- 
tion is evident from this plot. 



Figure 7. This the corresponding energy spec- 
trum of gas within the shocked cloud. Like the 
PDF, this spectrum too shows a distinct segreg- 
ation. The dotted lines show respective power- 
law fits; see text for description. 



Fuller (2008), Nutter & Ward-Thompson 2007, Motte et al. 1998). This proposition 
makes an investigation in to the origin of the CMF even more pertinent. Numerical 
simulations like those discussed by Vazquez-Semadeni (1994), and Padoan & Nord- 
lund (2002), among several other authors, and the PDF plotted in Figure 6 demon- 
strate the crucial role of interstellar shocks in generating a lognormal density PDF; 
although we note that the plot shown here is only semi-lognormal. Interestingly, a 
lognormal fit has also been attempted for the distribution of masses of dense (> lO'* 
cm"^) cores in nearby star-forming clouds (e.g. Enoch et al. 2008). While the ap- 
parent similarity of the dense-phase PDF derived here, with the CMF may perhaps be 
only a fortuitous coincidence, the likelihood of a causal relationship must be further 
investigated. Also, in the dense phase, as can be seen in Figure 6, this PDF peaks at 
~ 10^ cm which at 15 K, the precollision temperature within the cloud, corresponds 
to a characteristic mass, Mj ~ 3 M©; incidentally, using a gas-to-star conversion effi- 
ciency of 10% this characteristic mass corresponds to 0.3 Mq, the mass at which the 
stellar IMF develops a knee before turning over in to the brown-dwarf regime. This 
test simulation therefore seeks to re-emphasise the importance of interstellar shocks 
in forming putative star-forming pockets within molecular clouds, and may hold the 
key to explain the CMF. 

Power-spectrum The energy density, E(k), of a turbulent velocity field in the wavevector 
{k) space, the Fourier domain, is related to its Cartesian counterpart through a simple 
integral. 



The integral on the left-hand side extends over all possible wave-vectors, while the 
summation on the left-hand side runs over all SPH particles. The resulting power- 
spectrum for gas within the shocked cloud has been shown in Figure (7), which like 
the PDF, is also segregated. Although the spectrum for both, the dense as well as 
the rarefied phase of the gas, is a power-law it is, however, considerably steeper (oc 




(9) 
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kr'^-^) for the former, and Kolmogorov-like (oc ^ '■^) for the latter. The Kolmogorov- 
spectrum, as is well known, applies to an inviscid, incompressible fluid, though the 
gas here is viscous, but the rarefied phase appears to obey this approximation; it ob- 
viously breaks down for the dense phase. 

6. Conclusions 

1. I have argued in favour of my hypothesis that interstellar shocks could pos- 
sibly lead to formation of dense clumps, and contiguous filaments in molecular 
clouds. 

2. It has been demonstrated that propagation of a shock renders the density field 
unstable, and generates structure in it on a rather short timescale. Stars Ukely to 
form in these pockets, via energetic feedback, may further inject energy within 
the gas and quench star-formation in one region, only to trigger it elsewhere. 

Ancillary remarks This test case only examined the balance between self-gravity 
and thermal pressure as the magnetic field was not included, however, magnetohydro- 
dynamic simulations by for e.g. Padoan & Nordlund (2011) suggest a considerable 
modulation of the rate at which protosteUar objects form in turbulent gas. Other issues 
deserving a brief explanation include the eflFect of numerical resolution, and hydro- 
dynamic instabilities on the shocked cloud. The latter is essentially related to the 
numerical resolution, and more criticaUy, to the influence of SPH viscosity on dy- 
namically unstable fluid layers. First, a relatively poor resolution tends to suppress 
dynamical fragmentation. Second, the shearing interaction between the slab and the 
surface of the cloud is also likely to be unstable to hydrodynamic instabilities such 
as the Kelvin-Helmholtz instability. The demand on resolving this thin layer, how- 
ever, is conservatively large and not fulfilled here. It does not though, compromise 
the arguments presented in favour of the hypothesis examined here. 
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